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ABSTRACT 

We have reanalyzed the 900 ks Chandra X-ray spectrum of NGC 3783, findin g 
evidence on the asymmetry of the spectral absorption lines. The lines are fitted 
with a parametric expression that results from an analytical treatment of radia- 
tively driven winds. The line asymmetry distribution derived from the spectrum 
is consistent with a non-spherical outflow with a finite optical depth. Within this 
scenario, our model explains the observed correlations between the line velocity 
shifts and the ionization parameter and between the line velocity shift and the 
line asymmetry. The present results may provide a framework for detailed testing 
of models for the dynamic and physical properties of warm absorber in Seyfert 
galaxies. 

Subject headings: AGNs — absorption lines — ionized gas 


1. Introduction 

Since the launch of Chandra and XMM-Newton our view of active galaxies (AGNs) 
in X-rays has changed dramatically. It is now common to see rich and complex spectra in 
emission and absorption in the 0.5—10 keV band, indicative of absorption by partially ionized 
gas. The warm absorber hallmark, which appears in about 50% of AGN spectra, has been 
recognized since the Einstein observation of QSQ MR 2251+178 (Halpern 1984) and the 
subsequent detection of the O vm (871 eV) and O VII (739 eV) absorption edges by ROSAT 
and ASCA (Reynolds 1997; George et al. 1998; Nandra & Pounds 1992; Fabian et al. 1994). 
Then, the complexity of AGN spectra become more evident thanks to the improvement in 
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resolution and sensitivity of the instruments on board the new generation of X-ray telescopes; 
for example the HETGS (High-Energy Transmission Grating Spectrometer; Canizares et al. 
2000) with resolution of 0.012—0.023 A on board Chandra. This allows us to study in detail 
the physical structure and the dynamical conditions of the nuclear environment in Quasars 
and Seyfert galaxies. The best example of this so far is ngc 3783, which shows edges and 
absorption lines coming from H- and He- like ions of Si, S, Mg, Ne, O, N as well as L- and 
M-shell Fe transitions (Kaspi et al. 2002). 

One of the primary goals of all AGN studies is to gain an understanding of the gas 
flows in the vicinity of the central black hole. The X-ray absorption lines seen in NGC 3783 
are blueshifted with respect to the rest-frame of the host galaxy by moderate ~ 1000 km/s, 
velocity (Kaspi et al. 2001, 2002). UV observations with the Space Telescope Imaging Spec- 
trograph (STIS) on the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer 
(FUSE), have also revealed blueshifted absorption features coming from O VI, N v, C IV, 
N in and C III lines (Gabel et al. 2003a, b), providing a conclusive signature of outflow. This 
outflow character must be incorporated as an essential component in the unified picture of 
AGN, and is key to the understanding of the nature of these objects. 

Past work on the X-ray spectrum of NGC 3783 has yielded information about the warm 
absorber ionization, but has not provided detailed information or constraints on the gas 
dynamics. Fitting of photoionization models to the NGC 3783 Chandra spectrum requires at 
least two components with different ionization parameters in order to fit the line strengths 
(Kaspi et al. 2002; Krongold et al. 2003). There is evidence for kinematic differences between 
these components, but they are not apparent in all the lines. There may exist a dynamical 
relationship between the UV and X-ray absorbers. In NGC 3783 some of the velocity com- 
ponents seen in the UV are detectable in the X-ray lines, but not all (Crenshaw et al. 1999). 
There is no clear correlation between the velocity of the gas taken from X-ray measurements 
with the quantities associated with the ionization conditions in the gas such as the ionization 
potential of the ions (Kaspi et al. 2002) or ionization parameter. The picture that emerges is 
of an outflow with typical velocity 500-1000 km/s, but no clear evidence for the acceleration 
to these velocities, or about the geometry or driving mechanism for the flow. 

The 850 ksec Chandra observation of NGC 3783 is likely to remain among the most 
detailed and best-exposed spectrum of an AGN warm absorber for some time to come. Thus, 
some effort is justified in extracting the maximum information from it, and that is the main 
goal of the work we report in this paper. In doing so, we attempt to make use of a slightly 
different type of information contained in the spectrum than has been previously considered. 
Analysis efforts so far have focussed on understanding the profile of the absorption lines in 
this spectrum in terms of line centroids, strengths, and widths in both X-ray (e.g. Kaspi 
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et al. 2002; Krongold et al. 2003) and UV (e.g. Crenshaw et al. 1999; Gabel et al. 2003a). 
Kaspi et al. (2002) reported that most of the absorption lines in the X-ray spectrum of 
NGC 3783 show asymmetry, having more extended blue wings than red wings. However no 
tool has been developed to quantify this asymmetry line-by-line. One hypothesis which has 
been proposed to explain the asymmetries of the X-ray line profiles is that we axe observing 
two or more velocity components blended, with column densities smaller at larger velocity 
outflow (e.g. for the O VII lines). Support for this idea comes from the observation of multiple 
discrete components in the UV lines, some of which share outflow velocities with the X-ray 
lines. Recent observations show evidences of a correspondence between the UV velocity 
components and the X-rays (e.g., Krongold et al. 2003; Netzer et al. 2003; Krongold et al. 
2005). Although Krongold et al. (2003) and Netzer et al. (2003) among others, demonstrate 
that the inner shell transitions of Fe VII-XII seen in the X-ray region axe formed by conditions 
that favor the UV lines production, this is progressively less evident for X-ray lines coming 
from higher ionization states. Furthermore, the cross sections for the UV lines are greater 
than those for X-ray lines by factors > 10 2 3 , so that discrete components appearing in UV lines 
could result from very small fractional changes in the flow properties. Finally, it appeaxs that 
outflows from objects within our galaxy, such as hot stars, are characterised by essentially 
monotonic acceleration from subsonic speed to a speed which is characteristic of the escape 
velocity from the star. For these reasons, in this paper we will explore the hypothesis that 
the flow we observe in the X-ray band can be characterized as smooth and monotonically 
accelerated. This does not preclude the existence of discrete absorption components which 
affect the UV line profiles, but implicitly assumes that they are small perturbations on an 
essentially smooth flow. 

In this paper we present fits of the profiles of the strongest lines in the NGC 3783 
HETG spectrum using a parametrization technique similar to that used in profile fits in 
hot stars. This has the advantage that we can examine the implications of the fits for the 
global properties of the flow. There is clear evidence for asymmetry in some lines, and the 
line shapes contain information which can be interpreted in terms of the velocity law in the 
flow, the mass budget, and the ionization balance. We also find evidence for a systematic 
dependence of the profile shapes on the properties of the parent ion, and this has implications 
for the dynamics of the gas in the outflow from NCG 3783. 


2. Observations and data reduction 

The 900 ks spectrum of the Seyfert 1 Galaxy NGC 3783 used in this work is composed 
of six smaller observations, one of them reported by Kaspi et al. (2001) with a total 56 ks 
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of exposure time and the other five made in a large X-ray campaign (Kaspi et al. 2002) 
with ~ 170 ks each (see Table 1). We have analyzed all of these data uniformly using the 
most recent versions of the Chandra Interactive Analysis of Observations software (ciAO 3.0) 
and the matrix response calibrations (CALDV V.2.3) available to us to build new ancillary 
response files (arf files). We extracted the negative and positive first order spectra of the 
meg (Medium Energy Grating) and HEG (High Energy Grating). These have approximately 
0.023 A and 0.012 A spectral resolution, respectively. In the energy range 0.5—10 keV, the 
total numbers of counts in the first order is 583,196 for the MEG and 313,861 for the HEG. 
The S/N values at ~ 7 A are ~ 23 for the MEG and ~ 10 for the HEG in the combined 
900 ks spectra (Kaspi et al. 2002) using the default CIAO bins of 0.005 A and 0.0025 A 
respectively. We have checked all our measurements against the HEG data and find no 
significant diff erences. Moreover, we find the MEG data better because of the better signal 
to noise ratio in the spectral region of interest, and so in what follows we present results 
based on the analysis of the added negative and positive first-order MEG data. A detailed 
discussion of systematic uncertainties in wavelengths and flux measurements is presented by 
Kaspi et al. (2001), Kaspi et al. (2002) and references therein. Kaspi et al. (2002) reported 
a negligible background level with HEGTS/ACIS < 0.5 % of the signal. For this reason we 
did not subtract any background in this survey. In the analysis of the following sections we 
adopt single bin resolution elements, the default MEG size (0.005 A), giving us the maximum 
grating resolution. 


2.1. Fitting of the Continuum 

The continuum of NGC 3783 in the range 0.5 — 10 keV has been fitted by a power law 
with a photon index that spans from T ~ 1.5 to E ~ 1.83 (Blustin et al. 2002; De Rosa et al. 
2002), sometimes including a thermal component with kT ~ 0.1 keV (Krongold et al. 2003). 
Such fits depend on disentangling the true continuum from the effects of the absorption lines, 
which may include many lines which are unresolved and overlapping, and from the effects 
of possible emission. An accurate analysis of this type requires a model capable of fitting 
essentially all the absorption and emission features in the spectrum, including those which 
cannot be clearly resolved. Current spectral models have been very successful in fitting the 
majority of the strong lines (e.g. Krongold et al. 2003; Netzer et al. 2003), but currently 
available databases of atomic lines are insufficient to allow for highly accurate fits to the 
global spectrum (including weak features in both absorption and emission). For this reason, 
in our analysis we do not attempt to fit the global continuum. Instead, we work with distinct 
narrow spectral bands (~ 1 A) centered on the lines of interest. In each band we fit the 
spectrum to the line profile described below, plus a local continuum which is a power law 
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with T = 1.77 and normalization adjusted to each narrow band. 

In Figure 2 we show the continuum level linked together over the entire spectral range. 
In spite of the uncertainty associated with the continuum placement, we are confident that 
our choice of continuum does not introduce major errors into the fits of the absorption line 
profiles. This is because our continuum is a good fit to the local spectrum surrounding each 
of the strong lines, and because we have chosen our fitting intervals in such a way as to have 
the lines in the middle of each interval. This avoids dependence of the derived line properties 
on the continuum near the band boundaries, where there are unphysical discontinuities every 
~ 0.6 A, from the connections between some of the narrow bands. Nonetheless we take into 
account possible variations of up to 10% in the continuum level in our analysis of errors, as 
described in the next section. 


2.2. Line Identification and Atomic Data 

The spectrum of NGC 3783 contains a large number of absorption lines that are clearly 
apparent as local minims. Kaspi et al. (2002) identify 135 of them. In addition, the apparent 
bumpiness of the continuum suggests that there may be more weak or blended unidentified 
features. These complicate the fitting of the continuum which is reminiscent of the many 
emission lines of Fe II that affect the fitting of continuum spectra of AGN and some stars 
in the UV (Verner et al. 2003). In order to minimize the effects of blending and to allow 
for accurate determination of the outflow, we choose a subset of strong resonance lines with 
no obvious signs of blending (see Table 2), and whose wavelengths are accurately known 
(within ~ O.OlA) from experiment (e.g. Brown et al. 2002, and NIST 1 ). These lines were also 
identified by Kaspi et al. (2002) (see Table 3 of Kaspi et al. 2002). 


2.3. Line profiles 

In Figures 3-5 we show examples of several regions of the spectrum containing absorption 
lines. Several of these are clearly asymmetric. In the case of O VIII A 19 A there is also 
evidence for an emission component, reminiscent of P Cygni profiles. There is a few other 
lines in the spectrum which share this property, but in the great majority the emission is 
weak or absent. The goal of our analysis is to characterize and understand the dynamics 
of the gas flow from the shape of these lines, and in doing so we concentrate exclusively on 
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the absorption part. This is because the emission appears to be coming only from a small 
subset of the lines, and because the emission component in O VIII A19 A is narrower than 
the absorption component, so that the absorption at high Doppler velocities is safely out of 
the wings of the emission line. This approach has been shown to be fruitful in studies of 
the profiles in hot stars, for which the modeling of the emission component is less reliable 
than the study of the absorption (Lamers et al. 1987). Possible reasons for this include 
the contributions of thermal emission to the emission component, poorer understanding of 
atomic emission rates, and the fact that the emission depends more sensitively on the global 
properties of the flow geometry than does absorption. In AGNs it is likely that the flows are 
not spherically symmetric which will add further uncertainty to any theoretical prediction 
of the emission/absorption ratio. This provides further motivation to our decision to limit 
the study to solely absorption. 

Our procedure for analysis of the line profiles consists of fitting the flux in each narrow 
1 A band containing one of the lines of interest to a model of the form 

F(\) F :ont (A) x exp[-r(tw)J, (1) 


where A cont (A) is the continuum flux described in the previous subsection. 

Our wind model is based on the Castor & Lamers (1979) P Cygni profile atlas, built 
using the Sobolev approximation and parametrized by two functions, the velocity law v(r) 
and the optical depth as function of the velocity t(v). Then Lamers et al. (1987) developed 
the SEI method (Sobolev source function with Exact Integration) and applied to compute 
mass loss rates and terminal velocities in winds of Wolf-Rayet stars, being as accurate as 
the comoving-frame method of Mihalas et al. (1975) but much more numerically efficient 
to use interactively in the fitting of calculated profiles to the observations. In the Sobolev 
approximation the radial optical depth is (Castor et al. 1975) 


T S ob(v) 
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where / is the absorption oscillator strength, Ao (in cm) is the laboratory wavelength of the 
transition, rp, n u (in cm~ 3 ), gi, g u are the number densities and the statistical weights of 
the lower and upper levels of the transition respectively, and dv/dr is the velocity gradient 
in the wind. 


For resonance lines ( n u /g u <C ni/gi) one can approximate 
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The velocity law used throughout this work is based on the stellar wind theory of Castor 
et al. (1975), where the wind is assumed to be driven by the radiative force on spectral lines 
and its velocity varies with the distance as 


w(x) = w 0 + (1 - m 0 ) 



( 4 ) 


here w is the velocity normalized to the terminal velocity of the wind Vr&, w 0 is the velocity 
in the base of the wind and x is the distance normalized to the radius of the central core 
r/r o- The parameter (3 is the quantity governing the slope of the velocity with the distance 
and its ad hoc value depends on the type of radiative force acting on the wind. Analysis of 
hot stars suggests that 0.5 < (3 < 1, and we adopt the value /?=! (Lamers et al. 1987). 


^From this, the line optical depth can be parametrized as in Lamers et al. (1987): 


t{w) = (rtot//)(wM) Ql ( 1 “ ( w/wi) l,p ) a 2 (5) 


where w = v/voo = (c/v oo)(l — A/A 0 ). In this equation we use W\ = iq/^oo where Vy is the 
position of the blue edge of the line in velocity space (see, e.g., Lamers & Morton 1976). The 
total optical depth is 

T to t = ( ne 2 /m e c)f\ 0 Ni/v u (6) 

where TV* is the total column density of the absorbing ion, and 

I = wif y°‘ l {\-y 1,0 ) a2 dy (7) 

Jo 

is the normalization integral. 

It is important to note that although equation (5) was derived for radiative outflows, 
it is a convenient parametric formula for driving mechanisms as well, in which case the free 
parameters o l5 o 2 , and (3 may adopt slightly different physical meanings. 

In order to illustrate the role of each parameter on the line profile, we show in Figure 1 
the results of various models. In Figure 1(b), we can see the sensitivity of the profile when 
012 / 0:1 >> 1 (blue lines) for different values of o 2 . Figure 1(c) shows the case when a line is 
symmetric, i.e. (* 2 / 0:1 = 1, and Figure 1(d) shows lines with <* 2 / 0:1 << 1 (red lines). When 
<* 2 / 0:1 is large then the line has what we will call ‘blue skew’, that is the red edge is sharp 
and the blue wing is more gradual. In this case the slope of the blue wing is controlled by 
<*2 (with <*! fixed). When <* 2 / 1*1 is small then the line has what we will call ‘red skew’, that 
is the blue edge is sharp and the red wing is more gradual. In this case the slope of the 
red wing is controlled by ai (with <*2 fixed). This shows that the line profiles can be easily 
characterized by the value of 02 / 01 . 
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Fitting of lines in the spectrum of NGC 3783 was carried out with the Interactive Spectral 
Interpretation System package (ISIS ver. 1.1.8) using the xspec modules from the lheasoft 
libraries. In doing this we constructed a new XSPEC local-module named XSPCYG which is 
based on the computation of the optical depth as given by equation (5) (Lamers et al. 1987). 

Each line fit depends on four free parameters and one fixed; E ohs , v u on, a 2 and T tot . 
Sobs, is the energy of the red edge of the line (in.keV), which we take as the line transition 
energy in the rest frame of the galaxy. In Table 3 we present the results of the measurements. 
In the first and second col umn we have the ion and the rest wavelength of the line respectively. 
In the third column we give Vi, which is the position of the blue edge of each line (in km/s). 
The parameters describing the asymmetry of the line, on and a 2 , are in columns fourth and 
fifth. The integrated optical depth T tot (eq. [6]) is in the sixth column. Our determinations 
of the goodness of the fits were carried out computing a local xl- I n column seventh and 
eighth we give the number of bins included and our calculation of xl- After finding the 
best-fit parameters we compute the 90% of confidence level of all the parameters with the 
CONF command in ISIS. We included 10% on the variation of the continuum (through 
the normalization parameter) and V\ in the Monte Carlo searching, to take into account 
systematic uncertainties in the computation of the confidence intervals. We show in Table 3 
the limits associated to each fit parameter. These errors propagate as 

A (^) =— Aon-^A a 2 , (8) 

\a 2 J a 2 ol{ 


where Aon = an (max) — an (mm), Aa 2 = a 2 (max) — a; 2 (mm), and 
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where Am is a third of the bin size (0.005/(3A o )c km/s). 

In Figures 3, 4 and 5 we present the fits of three Lya: lines from H-like ions, which sample 
a range of ionization parameters, i.e. S XVI A4.729, Ne X A12.134, and O VIII A18.869. Two 
vertical dashed lines are drawn in these figures indicating the zero velocity point, i.e. the 
rest wavelength in the AGN frame, and where w/wi = 1/2 (i.e. where a perfectly symmetric 
line should have an = a 2 ). These figures are representative of the types of parameters 
which are found to fit the observed lines. They fall into two categories: (1) lines, like S xvi 
A4.73 A and Ne x A12.1 A (Figs. 3 and 4), which have a 2 — an, and (2) lines, like 
0 vm A19 A, which have a 2 an. The former are symmetric, or nearly so, and the latter 
are asymmetric, with a pronounced ‘blue skew’. In the case of O VIII A19 A, this line is 
rare in the sense it exhibit a clear emission component that accompanies the absorption, 
while for the vast majority of other lines the emission component, if present, is below the 
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noise level of the spectrum. Nonetheless, we use this line to try to estimate the errors 
introduced into the absorption line fits when emission component is ignored. To this end 
we include a Gaussian profile to mimic the emission component and obtain the best fit with 
the set of parameters: area=1.50 x lO'^tstxio - 6 photons/s/cm 2 , CENTER=19.147to;^4 A 
2 , SIGMA=9.26 x lO'^s'^xio-s A. worse case this emission has the effect of blue- 

shifting the minimum point of the absorption component by ~ 200 km/s (from ~ 780 km/s 
to ~ 980 km/s) and changes the a.\ parameter by up to ~ 20%, while a 2 that is mostly 
effected by the blue wing is affected much less. These a parameter change in the sense of 
adding a red skew to the line; on the other hand, most of the low ionization lines like O VIII 
exhibit blue skew. Thus, we conclude that even in the likelihood of emission components 
below the noise level in other lines in the spectrum, the overall results of the present analysis 
will be substantially unaffected. 

3. Analysis of line fits and the wind properties 

In the previous section we were able to fit the absorption lines in the spectrum based 
on model line profiles. In this section we check the distribution of the fits in the parameter 
space of Oj, a 2 > and w to explore systematic errors. Next, we study the relationship between 
the ionization parameters related to the absorption lines and the bulk flow velocity. This 
allows us to derive conclusions about the physical conditions of the gas flow in NGC 3783. 

3.1. Line asymmetry distribution 

The model that we have adopted for fitting the lines imposes an important requirement 
on the distribution of symmetry parameters a 1 and 02 on the velocity space. This condition 
is derived by finding the points of maximum absorption of every line, i.e. the minims of the 
absorption profiles, computed as 

where w is evaluated at the minimum of the profile in the observed spectra. It is important 
to point out that w is obtained directly from the spectra independently of line fits and the 
choice of the continuum levels. 



2 Is the observed wavelength. 
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From the line optical depth given by equation (5) we get 

f 1/13 _ a i/ a *P 

\wij l + oci/a 2 P' 


( 11 ) 


This provides a test on the consistency of the model and the line fits since it relates pa- 
rameters derived from the fits for each line independently, i.e. oq and a 2 , to parameters 
directly measured from the spectra without any fitting, i.e. w. In Figure 6 we plot the 
values of ( w/w\ ) 1 ^ with (3 = 1 vs. a\/a 2 for the lines considered in this work. For this 
plot, the error bars in w are taken as 1/3 of the bin size (0.005/Aoc km/s) and the errors in 
ct\/a 2 are derived from the fits shown in Table 3 and propagated as in equation (8). We find 
good agreement between the observed points and the behavior predicted by equation (11) 
indicated by the solid line in the plot. The computed correlation coefficient is r = 0.98 and 
considering weighted points x 2 — 34.63. Another important conclusion is that approximately 
90% of the points are in the range w/wi < 0.5 (i.e. ct\/ot 2 < 1). These blue lines are all 
formed in a region of the wind with velocities less than half the velocity reflected by the blue 
edge, Vy/2. These appear as the lines with the greatest ‘blue skew’, i.e. a 2 > oq. Only a few 
lines coming from similar ionization state (i.e. 12 < £ < 50) show the red skew character, 
e.g. Mg XI, Ne IX, and Fe XVII (see Table 3), where the uncertainty due to emission can be 
greater because of the ratio a.\/oc 2 > 1. 


3.2. The physical conditions in the warm absorber flow 


In order to consider the implications of the correlation shown in Figure 6 between 
(w/wi) 1 ^ and oq/oq, it is useful to discuss the relationship between the values of the param- 
eters which we derive in Table 3 and the parameters which describe the physical conditions 
in an AGN outflow. 

A simple and illustrative scenario for the origin of the warm absorber flow is a spherical 
wind centered on the source of continuum radiation where the continuum source provides 
both the illumination and the ionization of the gas in the flow. 

We describe the spatial position in the flow by the radius which we denote as r, and 
which we normalize to a scaled variable x = r/r 0 according to a fiducial position r 0 . The 
primary observationally derived quantity is the resonance line optical depth which can be 
written in the Sobolev approximation (eq.[2]) as (Castor et al. 1975) 


ire* 


TradC^) = /A 0 ni(r) 

me 


dv 

dr 


-l 
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( 12 ) 
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where / is the absorption oscillator strength, A 0 is the laboratory wavelength with n * the 
number density of the absorbing ions, and dv/dr is the velocity gradient in the wind. 

There are two functions responsible for the variation of the optical depth with velocity: 
the velocity gradient and the ion density. The former is computed from the velocity law 
given by equation (4). The ion density is given by 


rn{x) = A e x n H (x) x q it 


( 13 ) 


where A E = n E /nn is the abundance of the element with respect to H, = ni/n E is the 
ionization fraction, and Unix) is the gas number density 

n H (x ) = n 0 x~ 2 w ~ 1 (14) 

and 

= (^lL; r ° 2v ~) (15) 

where M is the mass loss rate and fj,mn is the average mass of the particles. At this 
point we need to adopt a model for the ionization balance. We assume that the flow is 
in photoionization equilibrium so that the ion fraction at every point depends only on the 
ionization parameter £ (the ratio of the ionizing flux F to the gas density n ff ). Our de fini tion 
of f is taken from Tarter et al. (1969) 


_ 47rF 

n H ’ 

and F = F 0 x~ 2 . By combining equations (16), and (14) 


( 16 ) 


£ = 


( 17 ) 


where £ 0 = 47rF 0 /n 0 . Finally we propose an expression for the ion fraction in the wind, based 
on the assumption that the ion fraction depends as a power law in the ionization parameter 

Qi ~ 9 o (£/ 0 7? , ( 18 ) 

where is the ionization parameter where q l is maximum, and qo ~ 1 is the maximum 
ionization fraction. Then the optical depth can be written 

r ra d(w ) = —f\ 0 no—q 0 A E w v+1/l3 ~ 2 . ( 19 ) 

me Voo 

Comparing with equation. (12) for the optical depth, this expression predicts a profile with 
0.2 = 0 and a\ = rj + 1/(3 — 2. If aq > 0 such a profile would had a very pronounced ‘red 
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skew’. This is in conflict with the results of our fits shown in Table 3 where we find ot^ > 
Qi and aq > 0. aq = 0 is ruled out for all the lines at the 2 a level, and for many at much 
greater level of confidence. An optical depth of this form also does not address the question 
of the mass budget, total optical depth, or the driving mechanism of the flow. 

The spherical outflow model likely provides an illustration of some of the physical effects 
important in a warm absorber flow, but clearly does not account accurately for the line profile 
shapes observed in NGC 3783. Other evidence for non-spherical flow geometry comes from 
unified scenarios (Antonucci & Miller 1985) in which the gas ~1 pc from the AGN center is 
concentrated in a torus which blocks a direct view of the center. If the warm absorber flow 
originates from the torus, as suggested by Krolik & Kriss (2001) and Chelouche & Netzer 
(2003) then it may have axial symmetry, but likely will not resemble a spherically symmetric 
flow close to its origin. Similar arguments apply if the flow originates from an accretion disk, 
as has been suggested for broad absorption line quasar outflows (e.g. Murray et al. 1995). In 
either case it is likely that the flow divergence does not obey r ~ 2 dilution due to collimation 
by a confining medium, distributed sources of gas in the flow, magnetic fields, or transition 
from a planar flow into a spherical flow. Also, the radiation field may not be solely due to a 
point source which is diluted by geometrically, owing to attenuation or to a source which is 
spatially extended on a length scale comparable to the flow size. 

In an attempt to take these effects into account we modify the model equations as 
follows. Firstly, we write the density as: 

uh { x ) = no£~ 2+K tw _1 , (20) 

where a positive value of k implies that the gas flow dilutes more slowly than in a free 
spherical expansion, i.e. that there are sources of gas embedded in the flow, or that the flow 
is confined. A negative value corresponds to sinks of gas in the flow, or expansion of an 
initially confined flow in a flaring geometry. Secondly, we allow the radiation flux to have a 
form 

F = F 0 x x~ 2 - p , (21) 

where p is an index to mimic a deviation of flux from the pure geometrical dilution case. 
This is expected if the gas has a significant optical depth (if p is positive) or if there are 
sources of radiation embedded in the flow (if p is negative). Here and in what follows, we 
set yd = 1 for simplicity. Then we find that 

f = £ 0 w x (1 - w ) r , (22) 

where T = p + k, and then, comparing the resulting expression for line optical depth with 
equation (12) we find 


a >2 = Tr) — k = pp + k(t] — 1) 


(23) 
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and 

cci = p — 1. (24) 

The quantities p and k provide insight into the conditions in the warm absorber flow. The 
parameter p depends on the physics of ionization, and so will vary from one ion to another, 
according to the differing conditions in the gas. p and k are universal quantities affecting all 
the lines with comparable weight. 

A more detailed discussion requires examination of ionization balance in order to es- 
timate the values for p. We expect the ion fraction <&(£) to be a function with a single 
maximum at some value of £ = £, and decreasing monotonically at larger and smaller values. 
7 ] describes the logarithmic slope of the <&(£) curve, and it is well known that the absolute 
value of this slope increases with |£ — £|. The ionization parameter in the flow is given by 
equation (22). If p = k — 0 then £ is an increasing function of velocity. Otherwise, £ is a 
decreasing function of position (corresponding to 1/(1 — w)) at large position. At w = 0.9 
£ = £ 0 /10 p+K . Thus, it is likely that p is not constant through the flow for most ions. 

We will nonetheless consider briefly here a highly simplified scenario in which we assume 
that p is constant for each ion throughout the flow. This is intended to illustrate qualitatively 
how systematic differences in line asymmetries can arise, not to test a specific dynamical 
model. Also for the following analysis we assume that the ion number density behaves as 
a power law with the ionization parameter with positive p. This is plausible because the 
maximum absorption depends on the product of with the gas density shifting this peak. 
As result we can make the analysis in the region where £ < £, (i.e p > 0). 

We assume here that the flow has a characteristic ionization parameter £o which is much 
greater than £ for all the ions in the flow. And we will assume that for each ion p depends 
only on the relative values of £ and £ 0 . Specifically, if p = |log(£) — log(£ 0 )| 1 ' 5 for all ions, and 
if log(£o)=4, p= 1 and k= 0.5 then we predict 0 < a\ < 4 and 1 < a 2 < 8. This is crudely 
consistent with the results in Table 3. Furthermore, the results of in Table 3 show an inverse 
correlation of both aq and a 2 with log(£). Figures 7 and 8 show an and a 2 plotted vs. log(£) 
as error bars, along with the best fit linear regression lines, plotted as long dashed curves. 
The values of £ were taken from the ionization balance curves of Kallman & Bautista (2001). 
The regression fits were calculated using the weighted least squares algorithm of Akritas 
& Bershady (1996), which takes into account the non-uniform errors in finding the best fit. 
The Pearson’s correlation coefficients are —0.49 and —0.55 for on and a 2 , respectively, versus 
log(£). Also shown in Figures 7 and 8, as the short dashed curves, are the values of oq and 
a 2 predicted by equations (24) and (23), with the formulation for p described above, p=l 
and /?=0.5. This shows that such an analytic model is an adequate description of the best 
fit regression line in both cases. What this shows is that the lines in NGC 3783 never have a 
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red skew; they either have a blue skew or are more nearly symmetric. The amount of blue 
skew is greater for lines with lower £, i.e. lines whose parent ion is most abundant at lower 
ionization par am eter. This can be explained, crudely, if the lines from low £ ions have a 
greater value of rj, i.e. they are more sensitive to changes in the ionization parameter of the 
flow, than lines from greater £ ions. 

Another way of representing the global properties of the outflow, is by examining the 
relation between ionization parameter at the absorption line minims vs. the line velocity. 
Prom our model such a relation is given by equation (22). In Figure 9 we plot the expected 
ionization parameter for maximum optical depth of each line vs. the measured reduced 
velocity at the line minimum, i.e. w. For this plot we have adopted a final wind velocity 
Voo = 1000 km/s. The figure also shows the curve given by equation (17), with (3 — 1 and F 
and £ 0 adjusted through y 2 fitting to the data. The values found for these two parameters 
are T = 5.18 ± 2.28 and £ 0 = 15410 ± 120. 

The results presented so far in this section can be summarized as follows: A simple 
spherical flow conserving mass predicts ionization parameter increasing monotonically with 
velocity. If the ion fraction is a monotonically increasing function of ionization parameter, or 
even if it decreases more slowly than 1/(3 — 2, then the line will have a net red skew which is 
not consistent with the fitted line profiles. Departures from this scenario allow for blue skew 
since they allow the optical depth to decrease with velocity. This can occur if the ionization 
parameter decreases at large distance or high velocity owing to additional sources of gas or 
collimation of the flow, or due to attenuation of the ionizing radiation. The observational 
result that the blue skew is more pronounced for lines with low values of £ suggests that 
the ionization fraction is an increasing function of ionization parameter, and that the rate 
of this increase is greater for ions with low £. That is, these ions are more sensitive to both 
the increase and decrease in £ as the wind accelerates. 


4. Summary and discussion 

We have reanalyzed the 900 ks CHANDRA spectrum of NGC 3783. The purpose of this 
work was to explore the kinematics of the warm absorber in absence of a linear correlation 
between the Doppler-shifts of the lines and the ionization energy of the absorbing ion, as 
noticed previously by Kaspi et al. (2002). We fit to model line profiles similar to the SEI 
profiles of Lamers et al. (1987), developed to model hot star lines. The line profiles predicted 
by this model are generally consistent with the observed data and allow us to characterize 
the main lines in the spectra, in terms of total optical depth, terminal velocity of the wind, 
and symmetry parameters. 
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Our model for the wind also allows us to explain the relationship in the spectr um 
between the ionization parameter at which the lines absorb and the velocity shifts of the 
lines. Moreover, a theoretically predicted curve (eq. [22]) can be adjusted to fit the observed 
data. In doing this we find a value for the exponent T ~ 5. This exponent results from the 
sum of two other parameters that represent a deviation of the radiation flux away from the 
simple geometrical dilution and deviation of the gas density from normal spherical decay 
under hydrostatic equilibrium. The large value of T suggests a combination of a decay of 
the radiation flux with distance from the source and departure from spherical dilution of the 
gas density. 

The model also predicts a relationship between the velocity line shift, measured at the 
line minimum, and the line symmetry parameters. The observed data agrees well with such 
a relationship. This is interesting because it involves the line shifts which are measured 
independently of the fits of the lines and the symmetry parameters for the lines, which 
depend on the wind properties as well as on the micro physics of the absorbing ions. The 
results of this work do not directly address the issue of the dynamics of the wind, its origin, 
or the details of ionization properties. Our results try to provide a framework for testing 
detailed dynamical models. 

The present interpretation about the asymmetry of the absorption lines as being due to a 
continuum wind flow is different from the interpretation of previous authors who propose that 
the X-rays might come from the blend of at least three UV components. But these two ideas 
may not be necessarily contradictory, because a continuum flow that yields a spectrum with 
line velocity spread between 500 and 2000 km/s could well be approximated by the emission 
of two of three discrete components, within the limitations (in resolution and signal-to-noise) 
of current spectra. However, the scenario of multiple components raises the question of how 
are these components related to each other and how the flow maintains the equilibrium. On 
the other hand, this scenario has the advantage of explaining the apparent kinematic relation 
between the UV and X-ray spectra. This point is currently under investigation within the 
wind flow model and will be the subject of a further publication. 


We thank the anonymous referee for constructive comments. The observation would 
not have been possible to analyze without the enormous effort of the Chandra team. 
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Table 1. Chandra Observation Log of NGC 3783. 


Sequence Number 

UT Start 

UT End 

Time (ks) a 

700045 

2000 Jan 20, 23.33 

2000 Jan 21, 16:20 

56.4 

700280 

2001 Feb 24, 18:45 

2001 Feb 26, 17:48 

165.7 

700281 

2001 Feb 27, 09:18 

2001 Mar 01, 09:10 

168.8 

700282 

2001 Mar 10, 00:31 

2001 Mar 11, 23:30 

165.5 

700283 

2001 Mar 31, 03:36 

2001 Apr 02, 02:48 

166.1 

700284 

2001 Jun 26, 09:57 

2001 Jun 28, 09:10 

166.2 


a Sum of good time interval. 


Table 2. Atomic properties of resonance lines in the X-ray spectrum of NGC 3783. 


Ion 

wavelength (A) 

lower level 

upper level 

f-value a 

S xvi 

4.729 

Is 2 S 

2p 2 P° 

4.129jE(— X) 

S XV 

5.039 

is 2 l S 

ls2p 3 P° 

7.610.E(-1) 

Si XIII 

5.681 

Is 2 *S 

ls3p 1 P° 

1.521jE(-l) 

Si Xiv 

6.182 

Is 2 S 

2p 2 P° 

4.136£(-1) 

Si xm 

6.648 

Is 2 *S 

ls2p 1 P° 

7.520£(— 1) 

Mg XII 

7.106 

Is 2 S 

3p 2 P° 

7.874£(-2) 

Mg xi 

7.473 

Is 2 X S 

ls4p !p° 

5.616£(-2) 

Fe xxill 

8.303 

2s 2 

2s4p 

1.404E(— 1) 

Mg xn 

8.421 

Is 2 S 

2p 2 P° 

4.142B(— 1) 

Mg xi 

9.169 

Is 2 4 S 

ls2p i P° 

7.383 15(—1) 

Nex 

9.708 

Is 2 S 

4p 2 P° 

2.896i5(— 2) 

Nex 

10.240 

Is 2 S 

3p 2 P° 

7.883£(-2) 

Ne IX 

11.547 

Is 2 *S 

ls3p 3 P° 

1.485£(— 1) 

Fe xxn 

11.770 

2s 2 2p 

2s 2 3d 

6.444E(-1) 

Nex 

12.134 

Is 2 S 

2p 2 P° 

4.147.E(-1) 

Fe xxi 

12.284 

2p 2 

2s 2 2p3d 

1.280£?(+0) 

Fe xviii 

14.373 

2s 2 2p 5 

2p 4 (3 P )3d 

3.151£(-1) 

Fe xviii 

14.534 

2s 2 2p 5 

2p 4 (3 p )3d 

2.051E(-1) 

O VIII 

14.832 

Is 2 S 

5p 2 P° 

1.393.E(-2) 

Fe xvn 

15.015 

2p 6 

2p 5 3d 

2.676iS(+0) 

O viii 

15.188 

Is 2 S 

4p 2 P° 

2.896S(— 2) 

O viii 

16.006 

Is 2 S 

3p 2 P° 

7.89115(— 2) 

O vn 

17.200 

Is 2 *S 

ls6p 1 P° 

1.509 JS7(— 2) 

O vii 

17.396 

Is 2 »S 

ls5p ip° 

2.679JS(— 2) 

O VII 

17.768 

Is 2 is 

ls4p !p° 

5.502E(-2) 

O vn 

18.627 

Is 2 IS 

ls3p 1 P° 

1.456S(-1) 

O viii 

18.969 

Is 2 S 

2p 2 P° 

4.151JS(— 1) 


a from xstar atomic data base 
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Table 3. Spectral line fits 3. 


Ion 

wavelength (A) 

v\ (km/s) 

0=1 

0=2 

Ttot 

# bins 

xi 

S XVI 

4.729 

2400+}™ 

i.02l$;S 

2 25"*” 4 ’ 27 

0.3981°;?®! 

19 

0.91 

S xv 

5.039 

2700+420 

°.90lJ;i° 

2 55+ 1 ’ 40 
z,oo _1.40 

0.3421°;?!° 

16 

1.30 

Si xiii 

5.681 

2100™ 

0 67+ 0 - 93 
U -°'-0.27 

1 43+ 2 72 
0. 14 

0.324l°J|? 

21 

1.47 

Si xiv 

6.182 

2100 +6O ° 

ziuu_ 1020 

0.681°°® 

z - u,3 -0.20 

0.8611°;?™ 

26 

1.42 

Si xiii 

6.648 

16201™ 

i-06lj;g 

1 54 4 " 0,84 
1 * O4 -0.84 

0.7391!™ 

20 

0.85 

Mg xn 

7.106 

16801™ 

i o 7 +1.92 
AaW -0.68 

2 71 4-1 ' 49 

0 39n"t° 150 

U.OZO_Q 143 

11 

1.39 

Mg xi 

7.473 

12601™ 

0.831J;™ 

n qc+°- 72 
u.oe_ 0 04 

n ooq+o 19 ^ 
u.zzw_ Q 103 

12 

1.33 

Fe xxiii 

8.303 

1500l®2 0 

0.05l°°« 

1 42+ 0,78 
■*•‘^—0.78 

n otro+0.317 
u.oo<5_ 0 054 

15 

3.26 

Mg xii 

8.421 

1740™ 

0 nn +0 - 70 
U - oU -0.20 

1 0Q 4 * 2,07 
i * uy -o.n 

0 847 4 ~ 0 ' 763 

u.cv±/_ 0 381 

14 

1.45 

Mg xi 

9.169 

16801®™ 

0 46 +0 ‘ 64 
— 0.22 

1 94+ 0 - 68 
i '^ 4 -0.68 

n 707 + 0.375 
* —0.288 

7 

3.24 

Ne x 

9.708 

13801™ 

H go+0.33 
O.-6O_ 0i0g 

0.271°;?® 

1 oi n^“0.653 
654 

16 

2.93 

Nex 

10.240 

37801™ 

1 QR 4 " 2 ' 77 
i.»o _ 0 96 

10.4511;™ 

D 97 Q + 0 - 208 

u.z^y_ 0 136 

27 

2.88 

Ne ix 

11.547 

13801™ 

2.651?;™ 

1 no+o.57 
i-l»_0.57 

0.9721°;!™ 

13 

1.92 

Fe xxn 

11.770 

19801™° 

o.45l°;°| 

1 7K+0-96 
'°-0.96 

0.6801°;?™ 

17 

0.85 

Nex 

12.134 

29401™° 

2.081J;?® 

z. 40 + 7-95 
°- 4S -3.02 

0 fi'SQ 4-0 - 296 
u.ooy _Q 2 96 

25 

1.43 

Fe xxi 

12.284 

1800l??° 

n 74+ 1 - 04 

U./4_o 37 

i.oell;!! 

1 91 O + 0.550 
j..ziz_ 0 680 

16 

1.35 

Fe xviii 

14.373 

2340lfS 0 

1 Q4 + 2 - 72 
x - y 4 -0.97 

5-3lll.ll 

0 430 4 ' 0 078 
u.^ou_ 0 387 

25 

1.41 

Fe xviii 

14.534 

11401™ 

O05+ 2 - 87 

2.43li;U 

1 oqo+ 0.826 

x.oyz _ 0 733 

14 

1.88 

O viii 

14.832 

27601™™ 

4.15l°|° 

14.4ll|;™ 

0.4501!;°?? 

29 

1.28 

Fe xvn 

15.015 

10201™ 

8.661°;^ 

1 70+O.95 

4.1731?;™! 

11 

1.05 

O viii 

15.188 

19801™ 

o.86iS;g 


0-5791!;??? 

22 

2.51 

O viii 

16.006 

21601°“ 

3-03111? 

8-19lo.ll 

0-6101!;®?? 

27 

1.45 

O vii 

17.200 

24001™ 

o 1 0 + 4.39 
616 - 1.56 

8.321UI 

0.4301!;??? 

28 

2.47 

O vii 

17.396 

27601™ 

o i c+3.03 

«-«il|;a 

0.6011!??? 

25 

3.14 

O vii 

17.768 

2280lgQ 30 

O 1 7+2.68 
z * i '-0.02 

4.80111? 

0.6781!;??? 

28 

3.86 

O vii 

18.627 

18001™ 

3-4ll?'i| 

6.3811;!? 

0 7n7 + 0 040 

13 

0.86 

O viii 

18.969 

37201^™ 

l-66lg;|| 

8-451?;!! 

n 700'*' 0,434 

u *' uu ~0.364 

44 

2.07 


a All uncertainties 90% of confidence computed with the command CONF in ISIS 
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Fig. 1. — Parameter roles: in A) we plotted velocity w vs 1/x (R/r) with the velocity index 
(3 = 1. In B) T tot jI = 1, f*i = 0.1 and the values of a 2 are labeled. C) shows symmetric 
profiles with T tot /I = 10 and a\ = a 2 =2 (black), 2.2 (red), 2.4 (green) and 2.6 (blue). D) 
Ttot/I = 1, <*2 = 0.1 and the values of a\ are marked. 
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Fig. 3. — Left panel: histogram data of S XVI A4.729 with its adjacencies in the observed 
wavelength space. The model and the continuum used are superimposed in red. Here a 
Gaussian has been included (see text for parameter values). Right panel: Velocity spectra 
with two vertical dashes lines at velocity zero and \/2w\. 
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Fig. 4. — Left panel: histogram data of Ne x A12.134 with its adjacencies in the observed 
wavelength space. The model and the continuum used are superimposed in red. Right panel: 
Velocity spectra with two vertical dashes lines at velocity zero and l/2w\. 
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Fig. 5. — Left panel: histogram data of Ovm A 18.969 with its adjacencies in the observed 
wavelength space. The model and the continuum used are superimposed in red. Right panel: 
Velocity spectra with two vertical dashes lines at velocity zero and \/2w\. 
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Fig. 8. — a 2 vs. log(£). Error bars show values from table 3, long dashed lines show best 
fit linear regression, short dashed line is analytic model. Note the differing vertical scale 
between from the previous figure. 
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Fig. 9. — Ionization parameter where the absorption is maximum £ vs the Doppler-shift 
velocity normalized to a terminal velocity of 1000 km/s w. The predicted model in red. 



